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We report on first principles calculations of the anisotropy of the intrinsic spin Hall conductivity 
(SHC) in nonmagnetic hep metals and in antiferromagnetic Cr. For most of the metals of this study 
we find large anisotropies. We derive the general relation between the SHC vector and the direction 
of spin polarization and discuss its consequences for hep metals. Especially, it is predicted that for 
systems where the SHC changes sign due to the anisotropy the spin Hall effect may be tuned such 
that the spin polarization is parallel either to the electric field or to the spin current. 
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Despite its relative smallness, the spin-orbit interac- 
tion (SOI) in solids gives rise to many phenomena of 
technological relevance and general scientific interest - 
well-known examples are magnetocrystalline anisotropy 
and anisotropic magnetoresistance. The anomalous part 
of the Hall effect [1| , which is observed in ferromagnets 
even in the absence of a magnetic induction field, re- 
sults from the spin-dependent transverse velocities which 
charge carriers acquire in the presence of a longitudi- 
nal electric field due to SOL In paramagnets the spin- 
dependent transverse velocities of spin-up and spin-down 
electrons are exactly opposite generating a transverse 
pure spin current, which is known as spin Hall effect 
(SHE) From the theoretical point of view SHE 

and anomalous Hall effect (AHE) are thus intimately re- 
lated and new insights into one of the two effects usually 
improve understanding of the other. 

While the SHE had been predicted theoretically al- 
ready in 1971 it was demonstrated experimentally for 
the first time in 2004 [6|. Since then the enthusiasm 
about the SHE has not abated. It has been studied ex- 
perimentally in semiconductors 0-0] and metals [l(| • As 
the SHE allows to access the spin degree of freedom of 
the electron without making use of magnetism it is be- 
lieved to play an important role in future generations of 
spintronic devices. 

It is well known both from theory and experiment (see 
Ref. [lH and references therein) that the anomalous Hall 
conductivity exhibits anisotropy, i.e., it is dependent on 
the orientation of magnetization M. Experimental evi- 
dence for the anisotropy of the SHE has been reported 
for AlGaAs quantum wells In contrast to the AHE 
the SHE is isotropic in cubic materials, i.e., in order to 
observe the anisotropy of the SHE non-cubic materials 
have to be considered (H)]. Besides bearing potential for 
applications the anisotropic Hall effects are also inter- 
esting from the point of view of information encoded in 
them about the Fermi surfaces and mean free paths of 
metals. Furthermore, it has been proposed [12j to ex- 
ploit the anisotropy in order to distinguish experimen- 
tally between inverse spin Hall effect and competing ef- 



fects caused by the magnetic field of the transport cur- 
rent. So far, ab initio calculations of the anisotropy of 
SHE have not been discussed in the literature. 

In the present work we undertake a detailed study of 
the anisotropy of the SHE in nonmagnetic hep metals 
and in antiferromagnetic Cr. In particular, the general 
expression for the orientational dependence of the SHC in 
hep and tetragonal metals is derived and the anisotropies 
are calculated from ab initio within the density functional 
theory. Generally, in antiferromagnets the SHE is ex- 
pected to allow the generation of pure spin currents like 
in paramagnets. While the AHE has recently been stud- 
ied in complex magnetic structures [l3j first principles 
calculations of the SHE have been limited to paramag- 
nets so far. Since the magnetic structure breaks the cu- 
bic symmetry antiferromagnets such as Cr always exhibit 
anisotropic SHE. 

The anisotropy of the AHE manifests itself in the de- 
pendency of the magnitude of the conductivity vector 
on the magnetization direction. The conductivity vec- 
tor <r AHE (M) relates the anomalous Hall current density 
J AHE to the electric field: 



I AHE 



= E x a 



AHE 



(M). 



(1) 



While the anomalous Hall current is always perpendicu- 
lar to the electric field, it is not necessarily perpendicular 
to the magnetization, since conductivity and magnetiza- 
tion vectors are not parallel in general. In the case of 
the SHE in paramagnets there is no magnetization vec- 
tor M to control, only the direction of the applied electric 
field can be varied. However, the spin polarization of the 
induced spin current depends on the direction in which 
the spin current is measured (see Fig. 1(a)). Hence, for a 
fixed electric field a given spin polarization S is measured 
only in a certain direction. Thus, in analogy to Eq. ([T]) 
we may write 



Q s = E x er(S), 



(2) 



where Q s is the spin current density and <x(S) is the SHC 
vector. If the magnitude of the SHC vector depends on 



the spin polarization direction S the SHE is said to be 
anisotropic. 

The spin current is characterized by velocity and spin 
polarization. Hence, the spin current density Q is a ten- 
sor in the 9-dimensional space R 3 ® R 3 spanned by the 
basis vectors ej C3 s s . For clarity we use the symbols 
s x , s y and s z to denote the unit vectors of spin polariza- 
tion while e x , e y and e z are the unit vectors of velocity. 
In addition to the SHC vector we define the tensor of 
SHC crfj, which has three indices: i denotes the direction 
of spin current, j the direction of applied external electric 
field, and s the direction of spin polarization of the spin 
current. The general expression for the linear response 
of the spin current density to an applied external electric 
field is given by 

Q = ^2<T s ij e i ®s s E :j . (3) 

ijs 

Comparing Eq. © and Eq. © we find that the SHC 
vector <x(S) and the SHC cr?- are related as follows: 

= 5 E e «*°«&. ( 4 ) 

ijs 

where Ok is the fc-th component of the conductivity vec- 
tor, S = (S x , S y , S Z ) T and ey s is the Levy-Civita symbol. 
Eqns. ([3]|4]) prove Eq. ([2]), which we conjectured above 
from analogy to Eq. ([T]). 

In cubic systems symmetry requires that erf^ = a xy eij S . 
Thus, the SHC may be expressed in terms of one material 
parameter, Eq. @ simplifies into Q s = uiLE x S, and 

the conductivity vector is <r(S) = cr* S. Since the mag- 
nitude of the conductivity vector, a xy , is independent of 
S, the SHE is isotropic in cubic systems. The relation- 
ship between the direction of spin current and the direc- 
tion of spin polarization in cubic systems is illustrated in 
Fig. 1(a). 

In contrast to cubic systems the SHE in hexagonal sys- 
tems is anisotropic. Consider the structure of the hep 
transition metals, which is illustrated in Fig. 1(b). If the 
electric field is applied along the x-direction, the magni- 
tude of the spin current in y-direction will generally differ 
from the one in z-direction since the x-axis exhibits only 
2-fold rotational symmetry. The spin current density in 
direction n = (0, cos 0, sin 0) T is 

n Q = — (a* y s z cos8 — s y sin 6)E X . (5) 

Note that according to Eq. ([3]) n • Q is a vector point- 
ing in the direction of spin polarization. We define the 
anisotropy of the SHE for spin polarization in the yz- 
plane as A zy — a xy — a\ x . For a general angle the 
components of the spin current with spin polarization 
parallel to n (Q\\) and spin polarization perpendicular to 



2 




FIG. 1: (a) Spin currents in cubic systems induced by an elec- 
tric field in z-direction. For electrons e moving in y-direction 
the spin polarization S points in ^-direction, while for elec- 
trons going in :r-direction the spin points in minus y-direction. 
(b) Hexagonal hep structure of the transition metal Ti. Due 
to the electric field in ^-direction a spin current flows along 
n = (0, cos 6, sin 9) T . In general n and S enclose an angle 
a 90°. 

n (Q±) are given by 

Qn = n • Q • n = -~A zy sm(26)E x , 

2 (6) 

Q± = Wlx+ &zy COS 2 6)E X . 

If A zy ^ 0, the spin polarization is perpendicular to n 
only if n is along the y or z-direction, otherwise spin po- 
larization and direction of spin current enclose the angle 
a = arctan(Qj^/Q||) ^ 90°, as shown in Fig. 1(b). It 
follows from Eq. © that is zero at the angle 



O = arccos J -a Z x/A zy (7) 

if a z xy and a v zx differ in sign. At this angle 9q the spin 
polarization and the spin current are collinear. This is 
an interesting constellation, which cannot occur in cubic 
systems. 

The case of spin current in x-direction and electric field 
E = (O,i?cos0,-Esin0) T in the yz-plane is simply re- 
lated to the previous one by a minus sign: The compo- 
nents of the spin current with spin polarization parallel 
and perpendicular to the electric field E are given by 
Oil = iA zy sin(20)£ and Q ± = -(a*. + A zy cos 2 B)E X , 
respectively. At the angle 0o, Eq. ([7]), the spin polariza- 
tion and the electric field are collinear. Thus, one can 
achieve collinearity of spin polarization and electric field, 
or collinearity of spin polarization and direction of spin 
current if a xy and a v zx differ in sign. 

If the electric field is applied along the z-axis, the same 
magnitude of the spin current will be measured in all 
directions perpendicular to the z-axis, since the z-axis 
exhibits 3-fold rotational symmetry. The spin current in 
direction n = (cos0, sin0, 0) T is in this case 

n • Q = {a x yz s x sin - a y zx s y cos 6)E Z . (8) 

Symmetry requires a v zx = tjy Z . Consequently, the magni- 
tude of the spin current is independent of and the spin 
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FIG. 2: (a) For the hep metals Sc, Ti, Zn, Y, Zr, Tc, Ru, Cd, La, Hf, Re and Os and for antiferromagnetic Cr the spin Hall 
conductivities Oy Z and a% y are shown as light (red) and dark (blue) bars, respectively, (b) Decomposition of the SHC of Sc into 
perpendicular and parallel components following Eq. (|6]). The angle a enclosed by the direction of the spin current and the 
direction of the spin polarization is also shown. At the angle #0=62.2° the component of the spin polarization perpendicular 
to the spin current vanishes and a=180°. 



polarization is perpendicular to both the electric field and ten as a Kubo formula: 



In the case of the hep structure the conductivity vector 
and the spin current density, Eq. ([2]), may be expressed 
in terms of the anisotropy as 



<r(S) = <7£ 2 S + (0,0,A*„&) T , 

= < 2 E x S + A zy S z (E y , -E X ,0) T . 



(9) 



Hence, only two parameters, a yz and A zy , suffice to de- 
scribe the SHE in hep nonmagnetic metals. This is a 
major difference to the anomalous Hall effect, where four 
parameters are needed in the phenomenological expan- 
sion of the conductivity of hep crystals up to third 
order in the direction cosines, because the band ener- 
gies depend on the direction of magnetization. Note that 
Eqns. ([5]|H]) apply also to tetragonal metals. 

For the sake of completeness we remark that the ana- 
log of Eq. and Eq. © runs J OHE = E x <r OHE (B) 
in the case of the low field ordinary Hall effect (OHE) 
in paramagnets, where B is the magnetic induction and 
<x OHE (B) = cr£ 2 B + (0,0, A zy B z ) T is the conductivity 
vector in hexagonal and tetragonal systems (cf Eq. ©). 



In general there is an intrinsic (independent of impu- 
rities) and an extrinsic (impurity-driven) contribution to 
the SHC (see Ref. Q and references therein for the origin 
of AHE, SHE is analogous). In the first principles cal- 
culations of the SHC presented below we consider only 
the intrinsic 0, [H 0, EH contribution to the SHC, which 
results from the virtual interband transitions in the pres- 
ence of an external electric field, and which may be writ- 



er, s ,, = eh 



d 3 k 



(27T) 



n;(k) = 23£ 



{^kn | Qj | V^km ) (4m | Vj | ^kn ) 



(10) 



where Vj is the j-component of the velocity operator, Q\ 
is the spatial i- and spin s-component of the spin current 
density operator, N occ is the number of occupied states, 
iV'kn) is the Bloch function of band n at fc-point k and 
s n is its energy eigenvalue. If only the spin conserving 
part of SOI is taken into account the spin current density 
operator may be written as Qf = h/2viT s . Here, t s is a 
Pauli matrix used to express the s-component of the spin 
operator. In order to treat the spin-nonconserving part 
of the SOI correctly we used the definition of the spin 
current density operator given in Ref. (lo |. 

To disentangle the intrinsic and extrinsic contributions 
to the SHC experimentally is still a challenge. The 
anisotropy of the extrinsic AHE is expected to be much 
smaller than the one of the intrinsic AHE 11] . Since 
AHE and SHE are analogous concerning their intrin- 
sic and extrinsic mechanisms, also the anisotropy of the 
extrinsic SHE is expected to be small. Thus, the di- 
rect comparison between the experimentally measured 
anisotropy of the SHC and the one calculated theoreti- 
cally based on Eq. (ITU)) allows to assess whether the first 
principles calculations predict the intrinsic contribution 
to the SHC quantitatively correctly for a given system. 
If quantitative agreement is found this provides a strong 
justification for the common procedure to attribute the 
difference between the experimentally measured SHC and 
Eq. (|TU1) to extrinsic effects, for the calculation of which 



ab initio methods have been developed recently 17-1 
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Our calculations of the intrinsic SHC, Eq. (ITOl) . for 
the hep metals Sc, Ti, Zn, Y, Zr, Tc, Ru, Cd, La, 
Hf, Re and Os and for antiferromagnetic Cr are based 
on the density functional theory and were performed 
with the full-potential linearized augmented-plane-wave 
(FLAPW) code FLEUR [20]. The generalized gradient 
approximation of the exchange correlation potential, a 
plane- wave cutoff of K max — 3.7 bohr , and the experi- 
mental lattice constants of the metals were chosen. In the 
case of Cr we neglected the spin-density wave and con- 
sidered the antiferromagnetic structure with two atoms 
in the unit cell and with the magnetic moments parallel 
and antiparallel to the z-axis. A dense fc-mesh is needed 
to perform the Brillouin-zone integration in Eq. (|10j) ac- 
curately. Consequently, we made use of Wannier inter- 
polation 21, [22| in order to reduce the computational 
cost. For this purpose we constructed a set of 36 max- 
imally localized FLAPW Wannier functions for each of 
the metals using the Wannier90 code (see Ref. [23| and 
references therein) and our interface [24] between FLEUR 
and Wannier90. 

The resulting SHCs are shown in Fig. 2(a). Except 
for Cd all metals studied in this work exhibit a large 
anisotropy of SHE, which we expect to be clearly visible 
in experiments. Of particular interest are the hep met- 
als Sc, Ti and Ru, where the sign of the conductivity 
changes as the spin polarization is rotated from the z- 
axis into the xy-plane. As discussed before, collinearity 
of the spin polarization and the electric field (of the spin 
polarization and the spin current) may be achieved if the 
electric field (the spin current) lies in the yz-plane at the 
angle 9q, Eq. ([7|), from the y-axis. To illustrate this we 
plot in Fig. 2(b) the angle a enclosed by the direction 
of the spin current and the direction of the spin polar- 
ization as well as the SHCs associated with Qn and Q± 
(see Eq. ([6])) as a function of the angle 9 for Sc. The 
critical angles at which the perpendicular component of 
the spin polarization vanishes are #o=62.2°, #o=32.1°, 
and 0o=19.1 o for Sc, Ti, and Ru, respectively. In the 
case of Cr the SHE is anisotropic as the cubic symme- 
try is broken by the staggered magnetization: If the spin 
polarization of the spin current is perpendicular to the 
staggered magnetization the SHC is larger by a factor of 
1.8 compared to the case of spin polarization parallel to 
the staggered magnetization. 

Generally, the integrand in Eq. (fT0"|) varies strongly as 
a function of k and the entire Brillouin zone has to be 
considered in the integration in order to reproduce the 
SHC quantitatively correctly. This makes it hardly pos- 
sible to interpret the spin Hall conductivity in terms of a 
small number of virtual interband transitions. Even the 
sign and order of magnitude of the SHC are difficult to 
predict based on simple arguments. Recently, the varia- 
tion of the sign of the Fermi-surface contribution to the 
SHC along the Ad and 5d transition metal series has been 
attributed to the variation of the sign of the spin orbit 



polarization on the Fermi surface [25|. However, we find 
that the spin orbit polarization does not change sign for 
Sc, Ti and Ru while the SHC changes sign as the spin 
polarization is rotated from the z-axis to the ?/-axis. 

In conclusion, we have investigated the dependence of 
the SHE on the directions of electric field and spin polar- 
ization. For the special cases of hexagonal and tetragonal 
metals we derived the general expression for the SHC vec- 
tor. We predict that in hep metals and antiferromagnetic 
Cr the SHE is strongly anisotropic. For Sr, Ti and Ru 
the anisotropy is particularly strong since the sign of the 
SHC depends on the orientation of spin polarization. In 
this case collinearity of spin polarization and electric field 
(or spin polarization and spin current) can be achieved 
for special directions of the electric field (or of the spin 
current). 
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[1 

[2 

[•3 
[1 

[5 

[6 
[T 

[8 
[9 

[10 

[11 

[12 
[13 

[11 
[15 

[16 
[17; 
[18 
[19 



N. Nagaosa, J. Sinova, S. Onoda, A. H. MacDonald, and 

N. P. Ong, Rev. Mod. Phys., 82, 1539 (2010). 

M. I. D'yakonov and V. I. Perel', JETP Lett., 13, 467 

(1971). 

J. E. Hirsch, Phys. Rev. Lett., 83, 1834 (1999). 

S. Murakami, N. Nagaosa, and S. C. Zhang, Science, 

301, 1348 (2003). 

J. Sinova, D. Culcer, Q. Nhi, N. A. Sinitsyn, T. Jung- 
wirth, and A. H. MacDonald, Phys. Rev. Lett., 92, 
126603 (2004). 

Y. K. Kato, R. C. Myers, A. C. Gossard, and D. D. 
Awschalom, Science, 306, 1910 (2004). 
V. Sih, R. C. Myers, Y. K. Kato, W. H. Lau, A. C. 
Gossard, and D. D. Awschalom, Nature Phys., 1, 31 

(2005) . 

J. Wunderlich, B. Kaestner, J. Sinova, and T. Jungwirth, 
Phys. Rev. Lett., 94, 047204 (2005). 
N. P. Stern, S. Ghosh, G. Xiang, M. Zhu, N. Samarth, 
and D. D. Awschalom, Phys. Rev. Lett., 97, 126603 

(2006) . 

S. O. Valenzuela and M. Tinkham, Nature, 442, 176 
(2006). 

E. Roman, Y. Mokrousov, and I. Souza, Phys. Rev. Lett., 
103, 097203 (2009). 

E. M. Chudnovsky, Phys. Rev. B, 80, 153105 (2009). 
S. D. Yi, S. Onoda, N. Nagaosa, and J. H. Han, Phys. 
Rev. B, 80, 054416 (2009). 

Y. Yao and Z. Fang, Phys. Rev. Lett., 95, 156601 (2005). 
G. Y. Guo, S. Murakami, T. W. Chen, and N. Nagaosa, 
Phys. Rev. Lett., 100, 96401 (2008). 
J. Shi, P. Zhang, D. Xiao, and Q. Niu, Phys. Rev. Lett., 
96, 76604 (2006). 

S. Lowitzer, D. Kodderitzsch, and H. Ebert, Phys. Rev. 
B, 82, 140402 (2010). 

M. Gradhand, D. V. Fedorov, P. Zahn, and I. Mertig, 
Phys. Rev. Lett., 104, 186403 (2010). 
M. Gradhand, D. V. Fedorov, P. Zahn, and I. Mertig, 
Phys. Rev. B, 81, 245109 (2010). 



5 



[20] See http://www.flapw.de 

[21] X. Wang, J. R. Yates, I. Souza, and D. Vanderbilt, Phys. 
Rev. B, 74, 195118 (2006). 

[22] J. R. Yates, X. Wang, D. Vanderbilt, and I. Souza, Phys. 
Rev. B, 75, 195121 (2007). 

[23] A. A. Mostofi, J. R. Yates, Y.-S. Lee, I. Souza, D. Van- 
derbilt, and N. Marzari, Computer Physics Communi- 



cations, 178, 685 (2008). 
[24] F. Freimuth, Y. Mokrousov, D. Wortmann, S. Heinze, 

and S. BKigel, Phys. Rev. B, 78, 035120 (2008). 
[25] H. Kontani, T. Tanaka, D. S. Hirashima, K. Yamada, 

and J. Inoue, Phys. Rev. Lett., 102, 016601 (2009). 



